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The visualization of curved sections of spacetime can be of considerable con- 
ceptual value. We discuss here the visualization of initial data for the prob- 
lem of the head-on collision of two black holes. The problem of constructing 



O ■ the embedding diagram is explicitly presented for the best studied initial 



data, the Misner geometry. We present a partial solution of the embedding 
diagrams and discuss issues related to completing the solution. 
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99; 1. Introduction 

> ■ 

^ . A benchmark for numerical relativity has been the computation of the gravi- 

tational radiation waveforms generated by the collision of two Schwarzschild 
throats starting from rest0. This description of the starting configuration 
is not complete. One can choose initial data in many ways to fit this ver- 
bal description. What remains to be specified can be said to correspond to 
the initial mutual distortion of the holes or the long-wavelength radiation 
present on the initial hypersurface. (Initial short wavelength radiation would 
presumably be "obvious.") When the initial separation of the throats is 
sufficiently small, a single nearly-spherical horizon surrounds both throats, 
and the exterior geometry can be thought of as a perturbation of a single 
Schwarzschild throat. Recent studies^] show that this viewpoint is reason- 
ably successful even when the horizon is highly distorted or is even split 
into two disjoint horizons. The extent to which the initial geometry mixes 
its "two throat" and "one throat" nature is the key to understanding this 
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perturbation approach which promises to be important in the future. One 
would like to get a feeling for how this changes with changing separation 
of the throats, with different choices of initial data for a given separation, 
etc. For these, and many other reasons, it would be very useful to have a 
direct way of visualizing initial data. Indeed, computer visualization would 
seem most appropriate for data that has been the starting point for so much 
intensive numerical computation. 

This would certainly seem possible in principle. The initial data under 
discussion are all momentarily stationary, so — in appropriate 4-dimensional 
coordinates — the initial time rate of change of the metric is zero; the initial 
data then consists only of the initial 3-geometry. Since that geometry is 
rotationally symmetric about the symmetry axis along which the holes will 
move, we can take a slice through that axis. The initial data is then fully 
specified by the spatial 2-geometry on that slice. A curved 2-geometry can, at 
least locally, be represented isometrically by a curved surface in flat Euclidean 
3-space. Due to the nature of the black hole initial data one supposes such a 
surface to have the general shape of a pair of trousers. Such pictures in fact 
are commonly drawn, but the actual surface, to our knowledge, has never 
been computed. 

In the present paper we discuss the generation of an embedding diagram 
for the most commonly used black hole initial data, the "Misner data" [0]. 
On rather general grounds we show that the generation of the embedding 
cannot proceed smoothly. We show explicitly that for the Misner geometry 
the breakdown takes the form of the mathematical equivalent of shock waves 
in the embedding surface. This breakdown occurs, however, only quite near 
the crotch of the trousers and does not stop us from computing most of the 
embedding surface. This "partial" embedding is sufficient for most aspects 
of visualization of the initial data. In a subsequent paper we will report 
on the extent to which it is possible to overcome the obstacles to a global 
embedding. 

The remainder of the paper is organized as follows: Section 2 intro- 
duces the mathematical preliminaries for the general problem of embedding 
a curved 2-dimensional geometry in Euclidean 3-space. Special emphasis is 
given to the case where the equations are everywhere hyperbolic since that 
turns out to be the case for the Misner geometry, which is described in Sec. 3. 
Numerical results are given in Sec. 4 along with a description of how the so- 
lution of the embedding breaks down at the formation of a "shock." In Sec. 5 
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we discuss whether these shocks are inevitable and, if they are, how the in- 
complete embedding diagrams can be of use. As a specific example, we show 
how they aid in the understanding of the range of validity of perturbation 
theory. 

2. The mathematics of embedding: 
The Darboux equation 

Consider a positive-definite, 2-dimensional geometry described in terms of 
local coordinates (x, y) and line element 

ds 2 = E dx 2 + 2F dx dy + G dy 2 . (1) 

We would like to realize this abstract 2-geometry as a curved surface in flat 
Euclidean 3-space, subject to the condition that the line element induced on 
the surface by the flat Euclidean geometry agrees with that of (|l|). We want, 
in other words, to find three functions (U, V, W) such that 

du 2 + dv 2 + dw 2 = E dx 2 + 2Fdxdy + G dy 2 (2) 

when u = U(x,y), v = V(x,y), w = W(x,y) are substituted into the left 
hand side of the above equation. This requirement leads to three conditions 

E = U 2 X + V 2 + W 2 X (3) 
F = U x U, y + V x V y + W, x W, y (4) 

G = K + K + K (5) 

where U }X means partial derivative of U with respect to x, etc. The functions 
(U, V, W) are called embedding funtions, and the mapping (x, y) — ► (U, V, W) 
defines a 2-dimensional surface in the (u,v,w) Euclidean 3-space. A set 
of functions (U, V, W) satisfying is referred to in the mathematical 

literature as a locally isometric embedding. Our goal is to find such an 
embedding for the 2-geometry corresponding to the Misner initial data. 

Rather than try to solve the above system of nonlinear, first-order partial 
differential equations (PDEs) for the unknowns (U, V, W), we will take a 
different approach originally due to Darboux. (See, for example, ||J].) This 
will lead to a single, nonlinear, second-order PDE for the embedding function 
W(x,y). Although the resulting equation is nonlinear in the second partial 
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derivatives of W, it is of a special type and can be reduced to a characteristic 
system of five quasilinear, first-order PDEs. It is this system of equations 
which we then try to solve numerically. Once we find W, the other embedding 
functions U and V are determined, in terms of W, by means of quadratures. 
The purpose of the remainder of this section (and the appendix at the end 
of the paper) is to make these statements more precise. 
Consider, instead of (|l|), the 2-dimensional line element 

du 2 + dv 2 = (E dx 2 + 2Fdxdy + G dy 2 ) - dw 2 . (6) 

Substituting w = W(x,y) into the right hand side of the above equation, we 
find 

du 2 + dv 2 = (E — W 2 X ) dx 2 + 2(F - W iX W, y ) dx dy + (G - W 2 y ) dy 2 . (7) 

Without loss of generality, we can assume that W jX and W >y vanish at the 
point about which we are trying to find the locally isometric embedding. 
(We can always perform a Euclidean transformation — i.e., a translation and 
rigid rotation — to guarantee that this is actually the case.) Since du 2 + dv 2 is 
positive-definite, and hence nondegenerate, in a neighborhood of that point, 
it follows that (TJ) must be nondegenerate for it to be be locally embeddable. 
The left hand side of (|7|) is also flat, and hence has vanishing Gaussian cur- 
vature. (There are no additional geometric constraints. In two dimensions, 
the Riemann curvature tensor is determined completely by the scalar curva- 
ture, and the Gaussian curvature equals one-half the scalar curvature.) If we 
express the Gaussian curvature of (|?|) in terms of the appropriate first and 
second partial derivatives of the components 

(E — W 2 X ) , (F — W tX W : y) , (G-W 2 y ) (8) 

we find 

A{rt - s 2 ) + Br + Cs + Vt + £ = (9) 

where A, B, C, V, £ are complicated functions of E, F, G and their first and 
second partial derivatives, and of the first partial derivatives p := W tX and 
q := W tV of W. (See the appendix for explicit expressions for A, B, ■ • • , £.) 
Here r, s, t are shorthand notations for the second partial derivatives W jXX , 
W }Xy , W m . All third derivatives of W cancel when calculating the Gaus- 
sian curvature of (0). A second-order PDE of the general form @, with 
A, B, ■ ■ ■ , £ independent of r, s, t, is said to be of the Monge- Ampere type. 
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It differs from the general, second-order PDE in that the only nonlinearity 
in r,s,t occurs in the combination (rt — s 2 ). This special feature results in 
a number of important simplifications below. In the context of the embed- 
ding problem — with A, B, ■ ■ • , £ depending on E, F, G as in the appendix — 
equation (||) is called the Darboux equation. 

If we write the general, nonlinear, second-order PDE for the unknown 
W(x, y) as 

Q(x,y,W,p,q,T,8,t) = Q (10) 
then the characteristic curves (x(X),y(X)) are those for which 

The discriminant 

A:=Q 2 s -AQ tr Q :t (12) 

determines whether real solutions exist for (x(X),y(X)). An equation is said 
to be hyperbolic for regions in which A > 0, so that two directions exist for 
characteristics at each point in the region. The equation is said to be elliptic 
for A < 0, and parabolic for A = 0. 

In the case of the Monge- Ampere equation 



Q,t 


= V + Ar 


(13) 


Q, s 


= C-lAs 


(14) 


Q ,r 


= B + At 


(15) 


= c 2 


- ABV + AAS . 


(16) 



yielding 

A 

The above expression for A does not contain r, s, or t. 

In the case of the Darboux equation, A has a remarkable geometric prop- 
erty: It can be written in terms of the Gaussian curvature K of the original 
line element in ([I]): 

A = -16K (EG - F 2 ) 3 (n 3 ) 2 . (17) 

(Here n 3 is the ^-component of the unit normal to the embedded surface. See 
the appendix for details.) For a positive-definite 2-geometry, EG — F 2 > 
and the sign of A is determined by the Gaussian curvature K. For K > the 



5 



Darboux equation is elliptic, so there are no characteristic curves. This means 
that features on one part of the embedding surface influence all other parts 
of the embedding. As we will show in the next section and in the appendix, 
the 2-geometry for the Misner initial data has everywhere-negative Gaussian 
curvature. In this case, the Darboux equation is hyperbolic, and features, 
such as the choice of boundary conditions for the embedding, propagate along 
characteristics and remain fairly localized. 

A geometric meaning for the characteristics of the Darboux equation 
emerges if we use (^3|)-([l5|) and the expressions for A, B, ■ ■ ■ , £ given in the 
appendix to rewrite ([□]) as 



-in'{EG - F 2 ) \k u f^V + 2K 12 K 22 

Here Kn, K12, K22 are the (x, y)-components of the extrinsic curvature tensor 
K for the embedding (see the appendix). The characteristic directions for 
the 2-geometry of ([!]) are thus seen to be the zero vectors (i.e., asymptotic 
directions) of K; they are the vectors v for which v ■ K • v = 0. Clearly, for 
K < there are two principal curvatures of opposite sign. The principal 
directions are orthogonal in the geometry of ([!]) and it is easy to see that 
the characteristic directions must be symmetrically arranged with respect to 
the principal directions. That is, a principal direction must bisect the angle 
between a pair of characteristics. 

Since we are guaranteed that there exist two characteristics through any 
point, we can use the characteristics themselves as coordinates. We label 





one family of characteristics (i.e., one family of solutions of (TT|)) with a and 
the other family with (3. The Darboux equation can then be reformulated 
in terms of five unknown functions (x, y, W, p, q) of the variables (a,/3). It 
is shown in the appendix that these are determined by the five quasilinear, 
first-order PDEs 

-(C-5)x >a -By >a -Aq, a = (18) 
^(C + 5)x > p-By iP -Aq,p = (19) 
^(C + 5)y, a -Vx, a -A P , a = (20) 
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~{C-8)y >f3 -Vx ! f 3 -Ap ! f 3 = (21) 
W, a -px, a -qy ja = (22) 

where 5 := v^A- (It should be noted that the general, nonlinear, second-order 
hyperbolic PDE, in characteristic form, requires eight equations for the eight 
unknowns (x, y, W,p, q, r, s, t). It is the special nature of the Monge- Ampere 
equation that results in a system of only five equations ||.) 

Of the three embedding functions (U, V, W), the Darboux approach sin- 
gles out one, W. It turns out, as shown in the appendix, that the remaining 
steps to complete the embedding are fairly straightforward. Once the so- 
lutions to fll8|) -(]22P are found, the remaining embedding functions U, V are 
determined in terms of W by means of quadratures 0. 

To close this section, we point out that the method of Darboux described 
above is only one of many approaches to solving the embedding problem. 
There is an extensive Russian mathematical literature on locally isometric 
embeddings for negative Gaussian curvature 2-geometries which describes 
the method of Riemann invariants. Interested readers should see the review 
article by Poznyak and Striking for more information. Also, a paper by 
Bernstein p[ describes an iterative numerical scheme for computing isometric 
embeddings. 



3. The Misner geometry 

The black hole initial data used for numerical relativity studies is the 
Misner 3-geometry|3[ given by 

rfs Misncr = ^ ^Misner [ V + drf + SU1 2 T] d(f) 2 ] (23) 

where 

OO 1 

V^Misncr = X/ / ' (^) 

n=-oo •v/COsh(yU + 2n/i ) — COS 7/ 



The geometry has a total ADM mass 



OO 1 

M « - ia £ s« (25) 
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and describes two throats located near // = ±/io- A measure of the separation 
of the throats is L, the proper distance from ji = —fi Q to \i — +/io, which 
can be shown to be 

/ 00 n \ 

L = 2a 1 + . (26) 

V ^sinh(n^o)y 

The constant a is a scaling factor that sets the size of both M to t and L, but 
not the ratio L/M tot , and does not affect the "shape" of the geometry. The 
only parameter affecting the shape is /zo, which is an increasing function of 
L/Mtot- For jjLo <C 1 there is a single, nearly spherical, initial horizon; for 
fiQ ^> 1 the geometry represents two widely separated throats, each with an 
initial horizon. The transition from a single horizon to a split horizon with 
two disjoint segments occurs at fi « 1.8. 

The Misner 3-geometry is rotationally symmetric — i.e., it is independent 
of the the azimuthal angle 0. We therefore lose no geometric information by 
taking a <fi = const slice of (f^|), thereby arriving at a 2- geometry 

and the possibility of an embedding in Euclidean 3-space. 

A numerical study of the Gaussian curvature K of this 2-geometry ([Z7j) 
shows that the Gaussian curvature is everywhere-negative. We have also 
been able to write the Gaussian curvature in a form in which it is manifest 
that it is everywhere-negative: 

K = ~ 2 a 2 ^Misner 

cosh (2p,o(m — n)) — 1 



E E 



m= -oon=-oo [cosh(/i + 2m/i ) — cosr]f /2 [cosh(/i + 2rijj, ) - cos?7] 3/2 

(See the appendix for details.) 

For our analysis, it is convenient to transform from the (//, n) coordinates 
to coordinates better suited to the description of the geometry at large dis- 
tances from the holes. To do this we transform from (//, rj) to (6,R), as if 
from bispherical to polar coordinates, according to 

= arctan (S:) (28) 




R = — : V sinh 2 fi + sin 2 77 . (29) 

cosh jl — COS Tj 
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The inverse transformation formulae are 



u = ±arccosh I R + a i ,^q\ 

/ (R 2 + a 2 ) 2 -(2aRcos9) 2 / 



rt = iarccos I . I . (31) 

'(R 2 -a 2 ) 2 + (2aRsm9) 2 



The signs are determined by requiring that fi > for —n/2 < 9 < +ir/2 and 
rj > for < 9 < Ti. In terms of the (9, R) coordinates, the 2-geometry (p7|) 
takes the form 

ds 2 = $ 4 (#, R) (R 2 d9 2 + dR 2 ) (32) 

with 

*=!+-£ 



R n+Q | sinh(n/io) | ^1 + (2a/ R) coth(nyU ) cos 9 + (a 2 / R?) coth 2 (n/i ) 

(33) 

If we identify (9,R) with (x,y), then the coefficients E,F,G of the line el- 
ement (ffl), which we take as the starting point for the embedding problem, 
are given by 

E = R 2 $\9,R) (34) 
F = (35) 
G = $ 4 (0,f2). (36) 

Since K < for the Misner 2-geometry, the corresponding Darboux equa- 
tion is hyperbolic and appropriate Cauchy data must be specified. We choose 
to specify "initial" data at large constant R. Here the 2-geometry ( |52] ) is ap- 
proximately that of a single, central Schwarzschild throat, and R plays the 
role of the Schwarzschild isotropic radius. R is related to the Schwarzschild 
curvature radius r curv via 

R = g r ™rv fl - M tot /r curv + yl - 2M tot /r rarv j (37) 

where M tot is the total ADM mass (^5|) of the single, central Schwarzschild 
throat. Equivalently, 

r curv = i? (1 + M tot /2R) 2 . (38) 
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For a single Schwarzschild throat of mass M to t, the embedding is given by 



W = 2M tot ^2R/M tot - ^/M tot /2Rj . (39) 

The initial values of W, p, q, r, s,t on the R = const surface are found from 
(p9|) and its derivatives with respect to 9 and R. 

Our numerical approach is a finite-difference solution to the characteristic 
equations in (|18|)-(|2~2"1). On the R = const initial value surface, grid points 
are chosen to be equally spaced in 6. A pair of (a, 0) characteristics is 
then started from each initial grid point. Equations fll8l) -(p2"D are used to 
propagate (x, y, W,p, q) forward along the characteristics. The two axes of 
bilateral symmetry of the Misner geometry are used to reduce the size of the 
numerical grid by a factor of four. The infinite sums were approximated by 
finite sums from —N sum to iV sum , with N sum large enough so that the omitted 
terms were negligible. 

4. Results: Characteristics and partial embed- 
dings for the Misner geometry 

Equations (|i8D-(^2[) were solved numerically, and from the results, solutions 
were constructed for U(a,/3), V(a,/3). With these solutions and the solution 
for W(a,(3), the two families of characteristics in the (u,v,w) space are 
known. The results are shown in Fig. 1 in the three cases Hq = 1,2,3 which 
correspond to L/M tot = 1.92,3.88,7.92. In these figures and in all those 
below, the value of the scaling constant a was chosen, for each value of //o, 
so that M to t = 1- The (u, v) coordinates therefore measure distances in units 
of M tot , so that in the figures we are always visually comparing spacetimes 
with the same total mass, but with throats at different separation. 

An immediately apparent, and crucial, feature of the figures is that the 
characteristic net does not cover the complete (u, v) interior to the R = 
const initial value surface. The numerical solution inevitably breaks down 
at the boundary of an oval region that includes the throats. This is not a 
numerical artifact; extensive numerical testing confirmed that the location 
of the breakdown was reasonably independent of numerical step size. 

A breakdown is, in fact, mathematically inevitable. This can most easily 
be seen by imagining the characteristics at the midpoint of the (u, v) plot. 
Since the underlying 2-geometry is symmetric about the rotation axis, the 
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embedding diagram must be reflection symmetric about the horizontal line 
through the midpoint. The symmetry axis, then, must be a principal di- 
rection of the extrinsic curvature. The zero directions of K, and hence the 
directions of the characteristics, must be symmetrically arranged about the 
principal direction. The nature of the characteristics must therefore be that 
shown in Fig. 2. But what direction along these characteristics is the "for- 
ward" direction of propagation of information from the boundary? A choice 
of direction along these characteristic segments violates the symmetry of the 
boundary. The direction of forward propagation must be ill-defined at the 
midpoint. An extension of this argument shows that this must be a problem 
not only at the midpoint but also, at least, on the symmetry axis joining the 
throats. The inescapable conclusion is that the characteristic net emanating 
from the R = const boundary cannot smoothly and continuously cover the 
entire region inside the boundary and outside the throats. 

This conclusion, in fact, does not depend on the symmetry of the Misner 
geometry, but only on its general topological character. We can shrink the 
R = const boundary, and the throats, so that the embedding region becomes 
a 2-sphere with three singular points, representing the boundary and each of 
the throats. Let us suppose that we can put a single family of characteristics 
on this 2-sphere, starting from the boundary and ending only at the throats. 
The tangent vectors to the characteristics (pointing in the direction of prop- 
agation of the solution) constitute a vector field having no singularities or 
zeroes except at the boundary and at the throats. At the boundary, the tan- 
gent vectors all point outward, while at the throats, they all point inward (if 
the specification of additional boundary conditions is to be advoided). The 
index of the vector field would then be +1 at each of these points, and 
the total index of the vector field would equal 3. But the Poincare index 
theorem requires that the index of a vector field on a manifold equal the 
Euler characteristic \ of the manifold. For a 2-sphere, the Euler characteris- 
tic x = 2, so the tangent field to the family of characteristics is impossible. 
A simple generalization of this argument shows that problems must develop 
in the embedding of negative Gaussian curvature initial data for any number 
of holes other than unity. 

It remains to determine the precise nature of the inescapable breakdown 
in the propagation of the characteristics. It cannot be due to a geomet- 
ric singularity. The intrinsic geometry of the embedded surface is fixed, of 
course, by the requirement that it be isometric to the nonsingular Misner 2- 
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geometry. We are carefully studying the extrinsic geometry of the embedded 
surface that we generate numerically, and we have found no evidence, so far, 
that any singularity is developing as the characteristics approach the point 
of breakdown. 

The actual nature of the breakdown is suggested by our argument above 
about the incompatible directions of characteristics. The breakdown occurs 
when characteristics of the same family cross each other, as in the classic 
example of shocks in gas flow. At such crossings there are three characteris- 
tic directions (two from the crossed characteristics, and one from the char- 
acteristic of the other family) and the propagation of the solution cannot 
proceed. Figure 3 shows a detail of the numerical results for the propagation 
of characteristics in the extreme case /xo = 6. To the left of the crossing at 
u ks 6.6, v ks 8.2 the results are valid; to the right, the results are meaningless. 

Although the development of these "embedding shocks" (i.e., character- 
istic crossings) block the computation of the embedding in an inner region 
of the geometry, we can still compute the outer region and obtain most of 
the visually useful information. Figure 4 shows the outline of the embedded 
surface (a view along the v axis) for the cases /xo — 1 and /i — 3. A clear 
geometric distinction is apparent. For the /xo = 3 case, the structure close 
in (in the region —5 < u < 5) is that of two throats. On a larger scale the 
geometry takes on the character of a single Schwarzschild throat. The shape 
of the /io = 1 case is quite different. Though it has the same "trousers" topol- 
ogy as that of /xo = 3, the steep sides of the embedding near u = ±2 mean 
that there is no region in which the individual throats look like more-or-less 
isolated holes. A 3-dimensional plot of the /x — 3 geometry is presented in 
Fig. 5. 

5. Discussion 

The development of "embedding shocks" is the most interesting feature of 
the problem of embedding the everywhere-negative Gaussian curvature Mis- 
ner 2-geometry. The primary question to be asked about them is whether it 
is inevitable that a "forbidden region" develop towards the crotch of the em- 
bedding, cut off from the outer region by the embedding shock. The topolog- 
ical arguments of the previous section clearly establish that something must 
go wrong with the characteristics, but the failure could be confined to the 
symmetry axis joining the throats, at which the characteristics meet with 
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incompatible directions of forward propagation. In effect, the "forbidden re- 
gion" could be degenerate, with zero area. This is not, of course, what our 
numerical results show, but the location of the embedding shock depends on 
the initial data. A correct small change in the initial data might shift the 
location of the embedding shock inward, reducing the area of the forbidden 
region; the precisely correct initial data might even reduce the forbidden re- 
gion to a line. We have numerically studied the effect on the shock location 
of the precise form of the initial data. The results strongly suggest that 
no improvement of the Cauchy data can significantly reduce the size of the 
forbidden region. [There is an additional point that must be added in con- 
nection with this. Part of the problem in the oval-shaped breakdown may 
be a technical difficulty that develops when the generated surface becomes 
vertical. This would be most expected near the v = plane at the outer 
edge of the throats. This problem, however, is causally disconnected from 
the breakdown near the u = plane, so the breakdown cannot wholly be a 
simple technical flaw in the approach. We are also studying this point, and 
developing a numerical technique for avoiding this problem, which may be 
inherent to the Darboux approach.] 

Suppose that somehow or another we were able to reduce the size of the 
forbidden region to that of zero area. What would be the implication of such 
a degenerate embedding shock? It would mean that the embedding diagram 
is uniquely determined by the initial data on the (approximately) circular 
outer boundary. That information would propagate inward unambiguously 
creating the embedding surface in its wake, until it hits the central degenerate 
shock. Whether the surface generated from one side and that generated 
from the other meet smoothly at this line is not guaranteed. The degenerate 
embedding shock would then be a sort of zero measure failure of the local 
isometric embedding to be global. This leads us to point out that theorems 
on the existence of global embeddings are notoriously nonexistent. There is 
no a priori reason to suspect that a global embedding of the Misner geometry 
either does or does not exist. 

If a global embedding does not exist, then an attempt to start at an outer 
circular boundary and propagate inward is doomed to incompleteness. We 
could, of course, start with inner data, and propagate the embedding diagram 
outward. But to discuss why we choose not to do this, it is useful at this point 
to review what we are trying to accomplish with the embedding. To help 
in the visualization of the physics, the embedding diagram must have the 
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character of a single throat at large radii, and in some sense must represent 
the idea of two holes at smaller radii. If a global embedding is impossible, 
then a diagram started from the central region would be guaranteed, at large 
radii, to look nothing like the embedding diagram of a single Schwarzschild 
throat. While there might be some uses to a diagram which represents the 
inner regions of the Misner geometry, at least for the visualization of initial 
data for black hole collisions, it would seem much more important to have 
the desired features at large radii, than in the inner region. 

In the discussion above, we have assumed that the forbidden region sur- 
rounded by the embedding shock can be reduced to a line, but this may very 
well not be the case. What if the embedding shocks are unavoidable fea- 
tures that separate the outer region from a nondegenerate inner region? The 
mathematical parallel between such embedding shocks and the more familiar 
shocks in gas flow are revealing. In the simplest gas shocks, called "kine- 
matic shocks," the equation of continuity, along with a relation between 



propagation velocity and density, leads to a crossing of characteristics jlO 



One then argues that the mathematical description is only an approximation 
to the true physics, and that the approximation breaks down at the shock 
formation. In principle, one can turn to a more complete description of the 
physics, such as the Navier-Stokes equation, to arrive at a mathematical for- 
mulation without singularities. The simpler approximate mathematics will 
be adequate on both sides of the shock; the more complete mathematical 
description will tell us how to "cross" it. That is, how to make the physically 
appropriate match of the two approximate solutions across the shock. 

This sort of strategy is inapplicable to embedding shocks. The Darboux 
equation is not an approximation; there is not a more complete geometric 
description that can tell us how to cross the shocks "correctly." But there is 
another view of shock- crossing conditions that may be of use. Without ref- 
erence to a more complete description of the physics, the crossing conditions 
for gas shocks can be inferred from conservations laws. In the simple case 
of kinematic shocks, the crossing conditions follow when mass-conservation 
at the shock is invoked. The resulting solution (without reference to more 
detailed physics) is then understood to be a "weak" solution — i.e., one which 
contains singularities, but for which the differential equation is satisfied when 
it is integrated over any region, including the region of the shock front [[Kj . 

What then is the geometric analogy in the case of embeddings? This 
question is presently being investigated. For now we can only point to 
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some interesting speculations. We note that for a sufficiently smooth sur- 
face, with extrinsic curvature everywhere defined and continuous, there can 
be no shocks, whether crossed correctly or not. It is interesting, however, 
to consider an embedding function W(u,v) in which there are special curves 
along which the second derivatives of W are bounded but discontinuous. 
Such a surface would be smooth (i.e., it would have no creases) but along the 
special curves, the extrinsic curvature would be discontinuous. At these spe- 
cial curves the directions of the characteristics (which are constrained by the 
extrinsic curvature) would be multiple-valued. One might imagine that such 
a surface would be a local isometry to the Misner geometry except along the 
special curves, and it is tempting to think that this, or some related geomet- 
ric phenomenon, would represent a weak solution to the embedding problem 
and the analogy of gas shocks. 

For the present, we point out that the embedding shocks do not preclude 
a partial solution of the embedding problem, which is potentially useful for 
visual and physical insight. An example is the question of the applicability 
of perturbation theory to black hole collisions. 

If the black holes are initially close together (i.e., if [xq is sufficiently small) 
the initial geometry outside the horizon is nearly spherical. The highly non- 
spherical geometry inside the horizon does not influence the time evolution 
of the spacetime ouside, so the evolution, and the generation of outgoing 
gravitational radiation, can be treated as a problem in perturbations of the 
spherical geometry of a single Schwarzschild throat. An important question 
is: How small must /i be for perturbation theory to be valid? As fi decreases 
below ~ 1.8 the topology of the horizon on the initial Cauchy hypersurface 
changes from that of two disjoint spheres (one around each throat) to that 
of a single 2-sphere surrounding both throats. Only for /^o significantly less 
than 1.8 will the horizon be reasonably spherical. Yet we know, from compar- 
ison with supercomputer evolution of the fully nonlinear equations 0, that 
perturbation theory works rather well for values of /j,q even somewhat larger 
than 1.8. The explanation would seem to lie in the the fact that nonsphericity 
just ouside the horizon disappears down the horizon. In the theory of pertur- 
bations of the Schwarzschild geometry [pTf , outgoing radiation appears to be 
generated at the peak of the "curvature potential," located at r curv ps 3M tot . 
This radius, roughly, separates the regions in which waves go outward to in- 
finity and go inward towards the hole. A reasonable criterion for applicability 
of perturbation theory would then seem to lie in the answer to the question: 
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How nonspherical is the geometry near r curv m 3M tot ? 

The difficulty in answering this question is that it requires a surface 
whose sphericity can be evaluated. The horizon is a natural choice but, 
of course, cannot be located at the peak of the curvature potential. In gen- 
eral, other choices are not gauge invariant, being dependent on coordinate 
choices. An exception is the coordinate-independent shape of the embedding 
diagram. We can then ask: How nonspherical is the embedding diagram near 
r curv ~ 3M to t? More specifically, we can evaluate the shape of the constant 
W contours corresponding to r curv ps 3M tot and see how they deviate from 
circles. 

In Fig. 6 we present such contours for /io = 1.5,2,2.5. The error that 
perturbation theory makes in estimating radiation energy in these three cases 
is 0, a factor of 2, and a factor of 10, respectively. For each value of /i , several 
contours are shown with the properties summarized in Table 1 below. For 
those contours with r := \Ju 2 + v 2 ~ r curv ps 3M tot there seems to be a clear 
correlation between the success of perturabtion theory and the spherical (that 
is, circular) shape of the contours. 



/i = 1-5 


/i = 2 


/i = 2.5 


'"max 


'"min 


Ar/r max a 


''max 


'"min 


Ar/r max 


'"max 


'"min 


Ar/r max 


2.50 


2.36 


.056 


2.77 


2.38 


.141 


3.05 


2.45 


.197 b 


2.94 


2.86 


.027 


3.12 


2.87 


.080 


3.56 


2.91 


.183 


3.60 


3.56 


.011 


3.72 


3.57 


.040 


4.02 


3.61 


.102 


4.48 


4.45 


.007 


4.56 


4.47 


.020 


4.75 


4.52 


.048 



a ^ r = r max — r m i n ^incomplete curve 



Table 1: Noncircularity of constant W contours 
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Appendix 

A.l Coefficients of the Darboux equation 

The coefficients A, B, ■ ■ • , £ of the Darboux equation are 

A = -4{EG-F 2 ) (40) 

B = +2p(2GF y -GG, x -FG iy ) + 2q(-2FF y + FG, x + EG, y ) (41) 

C = +4p(FG, x -GE, y )+4q(FE, y -EG, x ) (42) 

V = —2p(2FF x — GE X — FE jV ) — 2q(—2EF x + FE X + EE :V ) (43) 

£ = +(E-p 2 )(E, y G, y -2F, x G, y + G 2 x ) 

+(f- P q){F. ,.(:., - /•;./;,. - 2 /•;„/•;„ - 2 /•;,//.., + • !/•;.,. /•;,,) 

+(G - q 2 ){E, x G tX - 2E :X F jy + E 2 y ) 

+2((E - p 2 ){G - q 2 ) — (F — pq) 2 )(2F, xy - E m - G, xx ) . (44) 

Here W(x, y) is one of the three embedding functions (U, V, W), and p and 
q are shorthand notations for the partial derivatives W iX and W jV . The func- 
tions E, F, G are the components of the 2-geometry that is to be embedded. 



A. 2 Gaussian curvature 



The Gaussian curvature K of a 2-geometry having components E, F, G is 
given by 



K 



2(EG — F 2 )(—E jyy + 2F jXy — G jXX ) 

+ E ( G % ~ 2F ,xG,y + E^yG^y) 

+F(E tX G tV — E )V G )X + AF )X F )V — 2E jV F ty — 2F :X G :X ) 



+G{E 2 y - 2E tX F ty + E :X G :X )]/A(EG - F 2 ) 2 



(45) 



It is an intrinsic property of the 2-geometry, independent of any embedding. 
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A. 3 Unit normal 



The unit normal n to a surface depends explicitly on the embedding func- 
tions / := (U, V, W). Up to a normalization factor, it is given by the usual 
Euclidean 3-space vector product of f jX and f >y . When adjusted to be of unit 
length and reexpressed using the unit normal can be written as 

n = /'* X f ' y . (46) 



VEG - F 2 
The (u, v, w;)-components of n are 

i _ {V X W, V - V V W. 



n 



VEG 


— F 2 




WyU x ) 


VEG 


— F 2 




UyV x ) 



(47) 



n 2 = ^ ^§= "^ ZL (48) 

n 3 = ^^=J^. (I!)) 

VEG - F 2 v ' 

A. 4 Extrinsic Curvature 

The extrinsic curvature K is a tensor on the 2-geometry that describes the 
shape of the embedded surface and depends on the emebdding functions 
/ = (U,V,W). In terms of / and the unit normal n, its (x, y)-components 
can be written as 

Ku = f jXX ■ n (50) 
K 12 = f, xy -n (51) 

K 2 2 = f,yy ■ n . (52) 

Here • denotes the usual Euclidean 3-space dot product of two vectors. 

A. 5 The hyperbolic Monge- Ampere equation as a char- 
acteristic system 

We outline here how the hyperbolic Monge-Ampere equation can be writ- 
ten as a characteristic system of quasilinear, first-order PDEs for the un- 
knowns (x, y, W,p, q). More details can be found in the text by Courant and 
HilbertlBB. 
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We start with the Monge-Ampere equation 

Q:=A(rt-s 2 ) + Br + Cs + Dt + £ = 0. (53) 

This is a nonlinear, second-order PDE for a single unknown function W(x, y). 
The only requirement on the coefficients A, B, ■ ■ • , £ is that they be indepen- 
dent of the second partial derivatives r := W iXX , s := W jXy , t := W m . They 
can depend nonlinearly on x, y, and W, and on the first partial derivatives 
P '■= W )X and q := W tV . 

Assume that the equation is hyperbolic, so that 

A := Q 2 S - 4Q, r Q tt > . (54) 

In terms of the coefficients A, B, • • • , £, we have 

A = (C-2As) 2 - A(B + At)(V + Ar) 

= C 2 -ABV + AAS. (55) 



Here we used ( |53"D to eliminate r, s, t to obtain the second equality. The first 
equality is actually an important "identity" that we will use later on in our 
analysis. We choose to write this identity in the form 

C - 2As + 5 2{V + Ar) 
2{B + At) ~ C-2As-5 ^ ' 

where 

5 := VA = VC 2 - ABV + AAS . (57) 

Since equation (p3|) was assumed to be hyperbolic, 5 is real and is taken to 
be positive. 

We next write down the solutions to the characteristic equation 

In terms of the coefficients A, B, ■ ■ ■ , £ , the solutions are 

fdy\ C-2As±S 



dx) 2(B + At) 



(59) 



Since 5 > 0, there exist two distinct characteristic directions at each point. 
We can integrate along these directions to obtain the characteristic curves, 
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which we label by two parameters a and (3. If we choose a and (3 to be 
constant along the curves having the directions (dy/dx) + and (dy/dx)-, re- 
spectively, then 

( dy\ = V£ anc i I = ^ (60) 

\dx J + x t p \dx J x >a 

Here a; and ?/ are to be thought of as unknown functions of the characteristic 
coordinates (a,/?). The switch from (x,y) to (a,/3) as independent variables 
is allowed wherever the Jacobian of the transformation x^y^ — x^y )CC ^ 0. 
Using (|59|) and (|6~0"1), we obtain two equations 

{B + At)y, a -^(C-2As-6)x0 = (61) 
(B + At)y, p - -(C - 2As + 5)x, p = 0. (62) 



These two equations together with the identity ( p6|) give us two additional 
equations 

^(C-2^s + % ia -(P + „4r):r a = (63) 
l -{C-2As-5)y, p -{V + Ar)x tP = 0. (64) 

These are four first-order PDEs in the unknowns x, y. But we must still deal 
with the unknowns W,p,q,r,s,t. 

To this end, consider the first- and second-order strip equations. They 

are 

W >a = px tOC + qy ta (65) 
Wfi = px ; fs + qy i/3 (66) 



and 



P, a = rx ja + sy :a (67) 

P,p = rx >p + sy iP (68) 

q, a = sx ta + ty ta (69) 

Q,I3 = sx if3 + ty : p (70) 
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respectively. These equations are needed to guarantee that 



p = W, x q = W >y r = W, xx s = W, xy t = W m . (71) 

By substituting the second-order strip equations into flSID-flBip we can elim- 
inate the dependence on r, s, t, and write the results as 

By >a + Aq, a -^(C-5)x >a = (72) 

By 3 p + Aq t0 -^(C + 5)x^ = (73) 

~(C + 5)y, a -Vx, a -Ap >a = (74) 

\{C-5)y^-Vx tP -Ap,p = 0. (75) 

These are four quasilinear, first-order PDEs in the five unknowns (x, y, W, p, q). 
(Quasilinear in the sense that the equations depend linearly on the first par- 
tial derivatives of (x, y, W, p, q).) In order to complete the system, we need an 
additional independent equation. We choose this to be the first-order strip 
equation 

W t(x -px ta -qy ta = 0. (76) 



The above five equations are equations (^8|) -(p2|) in the main text. This is 
the desired result. 



A. 6 Determining U and V via quadratures 

After W has been found as a function of (x, y) we must still determine the 
remaining embedding functions U(x,y), V(x,y). We could in principle solve 
©"© directly, but there is a simpler approach, again due to Darboux, that 
we use and outline here. (For further discussion, see [§].) 

Consider a general, positive-definite, 2-dimensional line element 

ds 2 = f n dx 2 + 2/i2 dx dy + f 22 dy 2 (77) 

with metric components denoted by fy. Then the Gaussian curvature K of 
(|77[ ) can be written as 

JfK=(B tX -A, y ) (78) 
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where / := /n/22 — (/12) 2 is the determinant of the matrix of components 
fij and 



A 
B 



2%/7 
1 

V7 



In 
hi 

hi 



fll,x + fll,y — fl2,x 
fll,y — f22,x + fl2,y 



Using flT8| ) and Stoke's theorem, we can write 

/ JfKdxdy= I {Adx + B dy) 
Jn v Jan 



(79) 
(80) 

(81) 



where the integral on the left hand side of the above equation is taken over 
some simply-connected 2-dimensional region Q. The integral on the right 
hand side is around the closed 1-dimensional boundary dfl. 

If K = at every point (x, y), the right hand side of fl8~l|) vanishes for all 
closed curves dQ. This means that there exists a function d(x,y) such that 







A and d . 



B. 



(82) 



One can then integrate equations directly to obtain an expression for d 
in terms of the components fy. The solution is determined up to an overall 
additive constant. Moreover, using (|82"D and the definitions ([79|) and (|80f) of 
A and 5, one can show that 



/11 sin0 



fl2 n 1 \/7 Q 

: sin 17 H 7= cos » 



7. 



11 



Q 

H 7= COS 17 



vTii 
v 7 / 



sin0 



/11 v /11 

This means that there exist functions U(x,y), V(x,y) such that 
U, = 



V,. 



hi sin 

fl2 



Vhi 
+yAicosi? 
/12 



sin0 + cos ;i 
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cos + — ?== sin i) . 
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(83) 
(84) 

(85) 
(86) 
(87) 
(88) 
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These equations can also be integrated directly, allowing us to obtain expres- 
sions for U,V in terms of the components fy. The solutions for U,V are 
determined up to overall additive constants. 

The freedom in choosing the additive constants for U, V corresponds 
to the freedom of performing a 2-dimensional Euclidean motion of the plane. 
Such a motion preserves the form of the line element du 2 + dv 2 . The additive 
constant for •d corresponds to a rigid rotation of the plane about the origin. 
The additive constants for U, V correspond to a translation. 

Given the above results, we now specialize to the case of the flat 2- 
dimensional line element (0). This line element has components 

(89) 
(90) 
(91) 

where p = W )X and q = W tV . The solutions for U, V that we find by 
integrating and (|85|)-(|38"|) (using the above /y) complete the embedding. 



/ll 


= E 


~p 2 


/l2 


= F 


-pq 


/22 


= G 


-q 2 



A. 7 Gaussian curvature of the Misner 2-geometry 

We start with the Misner 2-geometry 

ds 2 = aVLner + V] (92) 

written in terms of the coordinates To simplify the notation in what 

follows, we will drop the "Misner" subscript from y3 Misner and define 

[n] := cosh(/i + 2nfi ) — cost]. (93) 

Then 

oo 

V := ^Misner = H ■ ( 94 ) 



Our goal is to calculate the Gaussian curvature K of (P2|) , and to write 
it in a form which is manifestly everywhere- negative. If we take (/i, rj) as our 
coordinates (x, y) then 

E = G = a 2 <p 4 and F = (95) 
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are the components of the Misner 2-geometry. For E = G, F = 0, the 
expression (|45f) for the Gaussian curvature simplifies considerably: 



K = --E- 1 [(\nE), xx + (\nE) 



yy - 



(96) 



This is just minus one-half times the covariant Laplacian of \nE. For E 
a 2 ip 4 , we have 



K = -2a- 2 V - i [(ln<p)^ + (ln ¥ >) 



vv . 



(97) 



This is what we must evaluate. 

We begin by calculating the first and second partial derivatives of cp. They 

are 



OO J 

£ --M- 3/2 sinh(/i + 2n/i 



n=— oo 

OO 



1 



E -JN 3/2 sin 7] 

n=— oo 



1 



Yj ( ~ 3/2 cosh(^ + 2n/i C 



4 l 



n 



~ 5/2 sinh 2 (/i + 2n/i ) 



7/7/ 



E 



- [n] 3 ^ 2 cos 77 + ^ [n] 5 ^ 2 sin 2 77 



(98) 
(99) 



(100) 
(101) 



Since 

(In <p) iMl + (In (p) m = (f~ 2 [(f (f^ + <p ip m - {(f^) 2 - (tp iV ) 2 
it is also convenient to evaluate 



(102) 



V 5 V,uu + V V,: 



00 00 



1 ^ ~ cosh(/i + 2m/i ) cosh(/i + 2n/n ) - cos 2 1] 
4 „^rL [m] 3/2 [n] 3 / 2 1 j 



m=— 00 n=— oo 
\2 1 /,„ 



(^) + W 

1 ~ ~ sinh(/i + 2m/i ) sinh(^ + 2n/x ) + sin 2 ?] 
— 7 r™is/2 r„ia/2 • ( 1U4 ) 



m=— 00 n=— 00 



[ m ]3/2 [ n ]3/2 
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Combining ([I02D , (|T03|) , and (|TD|) we see that 
(ln^), w + (lnp) 



oo oo 



- -2 cosh(2/x (m - w)) - 1 

m=— oo n=— oo L J L J 



From this we get the final result 

K 1 -2 -6 v- v- cosh (2 Mo (m -,;))-! 

2 * m E„S ■ (1 ° 6) 

This is the manisfestly negative-definite expression for the Gaussian curva- 
ture given in the main text. 

One may rightfully worry about those places where 

[n] := cosh(/i + 2n/i ) — cost] = . (107) 

When this happens, (p = and K — > 0. This occurs only if 

77 = and fjt + 2nfi = (108) 

where n is an integer. For — /j < /j < /i , we see that the second equality 
holds only if /x = and n — 0. Thus, [n] = if and only if (/!,??) = (0,0). 
This is the "point" at spatial infinity, and there we expect the 2-geometry 
corresponding to the Misner data to be flat. 
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Figure Captions 

Fig. 1: The net of characteristics propagating inward from the R = const initial 
value surface for /i = 1, 2, 3. 

Fig. 2: The principal directions of the extrinsic curvature and the characteristic 
directions near the center. 

Fig. 3: The computed crossing of two characteristics of the same family for 
Ho = 6. 

Fig. 4: Shapes of 3-dimensional embedding diagrams when viewed along the v 
axis for fi = 1, 3. 

Fig. 5: Perspective view of the (incomplete) embedding diagram for /i = 3. 
Fig. 6: Constant W contours of the embedding diagrams for /i = 1.5, 2, 2.5. 
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